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Are the missing X-ray breaks in Gamma-ray Burst 
afterglow light curves merely hidden? 
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ABSTRACT 

Gamma-ray burst afterglow observations in the Swift era have a perceived lack of 
achromatic jet breaks compared to the BeppoSAX or pie-Swift era. Specifically, rel- 
atively few breaks, consistent with jet breaks, are observed in the X-ray light curves 
of these bursts. If these breaks are truly missing, it has serious consequences for the 
interpretation of GRB jet collimation and energy requirements, and the use of GRBs 
as cosmological tools. Here we address the issue of X-ray breaks that are possibly 
'hidden' and hence the light curves are misinterpreted as being single power laws. 
We do so by synthesising XRT light curves and fitting both single and broken power 
laws, and comparing the relative goodness of each fit via Monte Carlo analysis. Even 
with the well sampled light curves of the Swift era, these breaks may be left misiden- 
tified, hence caution is required when making definite statements on the absence of 
achromatic breaks. 
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1 INTRODUCTION 

The afterglow emission of Gamma-Ray Bursts (GRBs) 
is well described by the blast wave, or fi reball, model 
|Rees fc Meszarodll992l ; iMeszaros et al.lll99St ). This model 
details the temporal and spectral behaviour of the emission 
that is created by external shocks when a collimated ultra- 
relativistic jet ploughs into the circumburst medium, driv- 
ing a blast wave ahead of it. The level of collimation, or jet 
opening angle, has important implications for the energet- 
ics of the underlying physical process and the possible use of 
GRBs as standard candles. The signature of the collimation, 
according to simple analytical models, is an achromatic tem- 
poral steepening or 'jet break' at ~ 1 day in an otherwise 
decaying, power law light curve; from th e time of this b reak, 
the jet opening angle can be estimated (Rhoads 1997). 

Since the launch of the Swift satellite ( Gehrels et al.l 

120041 '). this standard picture of afterglows has been called 
into question by the lack of observed ach romatic temporal 
breaks, up to weeks in a few b ursts (e.g.. IPanaitescu et al.l 
l200d :lB urrows fc Racusmll2007l) . In some afterglows, a break 
is unobserved in both the X-ray and optical light curves, 
while in other bur sts a break is obse rved in one regime but 
not the other (e.g-. lLiang et al]|2007l ). In the BeppoSAX era, 
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most well sampled light curves were in the optical regime, 
while in the Swift era most well sampled light curves are in 
the X-ray regime. Our expectations of the observable signa- 
ture of a jet break are hence based on the breaks observed 
pie-Swift, predominately by optical telescope s, and th e mod- 
els w hich explained t hem, notably those of iRhoadsl (|l997l . 
Il999h and ISari etafl l| 19991 ). It is not clear that the breaks 
will be identical in both the X-ray and optical regimes. In 
the cases of GRB 990510 and GRB 060206, both have clear 
brea ks in the optical but only marginal bre aks in the X- 
ray (|Kuulkers et all l2000l : ICurran et ai1l2007l respectively). 
Regardless of this, GRB 990510 is taken as a prototypical 
achromatic jet break, while the achromatic nature of the 
GRB 060206 break is only evident when supported by the 
broad-band spectral indices. 



In this paper we address the issue of X-ray breaks which 
are possibly 'hidden' and hence the light curve misinter- 
preted as being a single power law. We do so by synthesis- 
ing X-ray light curves and fitting both single and broken 
power laws, and comparing goodness of each fit via the F- 
test. In 3D we introduce our method while in [J3]we present 
the results of our Monte Carlo analysis. In Sjjjwe discuss the 
implications of these results in the overall context of GRB 
observations and we summarise our findings in ij5] 
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Table 1. The values of tem poral indices, a, for the gi ven values 
of p in the given cases, from lZhang fc Meszaroa (|2004h . 



p 




1.8 


2.4 


3.0 


vx > Vc,m 




0.96 


1.30 


1.75 


Vm < fx < V c 


(Wind) 


1.23 


1.55 


2.00 


I'm < fX < 


(ISM) 


0.71 


1.05 


1.50 


post-break 




1.95 


2.40 


3.00 



2 METHOD 
2.1 Models 

In accordance with the fireball model, we ascribe the be- 
haviour of the X-ray light curve to be a single power law 
decay where the flux goes as: 

F v {t) <xt~ a , 

with a temporal decay index, a. Alternatively, if there is a 
break we use a smoothly broken power law decay with a 
break at tb- 



F u {t) = F„{t b ) 



aiS 



+ 



-1/S 



where S is the smoothing factor (a larger value correspond- 
ing to a sharper break) and the light curve goes as F v oct~ ai 
and F v oc t~ a2 pre- and post-break respectively (ai < o^)- 
We have simulated both single and broken power law 
data for a limited, but varied number of realistic parame- 
ter sets. The flux offsets are set at a reasonably high rate 
of O.lcts/s at 1 day to give well sampled light curves. We 
take break times at 0.35, 1 and 2 days, and input smooth- 
ing parameters, S = 1,2,5. The temporal decay indices 
are chosen as a function of electron energy distribution in- 
dex, p, assuming i) X-ray frequency, vx > Vo,m, ii) v m < 
vx < Vc in a wind-like environment and iii) f m < v x < Vc 
in an ISM- like (constant density) environment, where v c 
and Vm are the synchrotron c ooling and peak frequencies 
(e.g. IZhang fc Meszarosl l2004t) . We adopt three values of 
p = 1.8, 2.4, 3.0 covering the observed range of p values (e.g. 
IPanaitescu fc Kumarll200ll ; IStarling et al.ll2008l ). which im- 
ply the temporal indices, a as shown in Table[T] 



2.2 Synthesis of data 

In sy nthesising Swift X-Ray Telescope (XRT; lBurrows et al.l 
2005) li ght curves we wanted them to be comparable to t hose 
produced by the online repository (|Evans et al.ll2007l ), St..) 
we tried to satisfy the criteria that were applied to their 
production process. It is therefore instructive to describe 
briefly their method, and the differences to ours. Firstly, 
lEvans et al.l present a reduction method dealing with ob- 
served source and background counts whereas ours deals 
directly wit h the rate predicted by an underlying model. 
I Evans et al.l demand that each bin (i.e., data point) must 
have a minimum time span and number of counts from the 
source region. The minimum number of counts per bin is 
dependent on which mode the XRT is observin g in (i.e., 
Windowed Timing/Photon Counting (WT/PC): Hill et all 
2004, 2005) and the measured counts of the source region. 



As we will be dealing with late time data (> 10 4 s) which is 
predominately measured in the PC mode, and has low count 
rates (< 1 cts/s), we use a constant value for the counts per 
bin. We do this by setting the bin interval so that the ap- 
proximate numerical integration of the rate, i.e., the counts, 
is constant. Also at these late times, it is unlikely that this 
minimum count will be reached in less than the minimum 
time span so the time span criterion may be ignored. 

The main conditions of our synthesised light curves are 
as follows: 

• Constant counts, and la fractional error of 0.25, per 
data point 

• 94 minute orbits (47min on/off) 

• Fractional exposure of 0.1 after 1 day 

• Rate cut off at 5 x 10 -4 cts/s 

• Time spans 10 4 to 2.6 x 10 6 s (3 hours to 30 days) 

Our constant co unts per data point is not comparable to 
the constant used by I Evans et al.l . though we have chosen a 
value so that the number of data points over our light curve 
is comparable to XRT light curves from the repository. The 
error value of an individual data point is taken as a constant 
(since the counts per bin are constant) fractional value of 
0.25, to agree with the values o btained from a number of 
well sampled XRT light curves in lEvans et all 

As the Swift satellite is in a low Earth orbit with a pe- 
riod of 94 minutes it suffers from 50 per cent time off target 
per orbit. We take this into account by not synthesizing data 
for these off-target periods, assuming that we start our light 
curve at the start of the on-target time. Even after this or- 
bital time loss, Swift does not dedicate 100 per cent of the 
time that the target is visible to the target itself. This is due 
to other observing commitments and constraints. In fact the 
fractional exposure or actual fraction of time on target after 
th e orbital con straint, is 0.1 after ~ lday for most bursts 
in lEvans et all leading to ~ 5 minutes of exposure time per 
orbit. 

As we are only interested in the power-law decay or af- 
terglow phase of GRBs, we start our synthesis at 1 x 10 4 s. 
This is a time after which many real XRT light curves 
have data unaffected by early flaring or prompt emission 
|Chincarini et all2007l ) . XRT can sometimes continue to ob- 
serve for up to a month (~ 2.6 x 10 6 s) so we synthesise data 
out to that same time, unless it drops below a threshold of 
5 x 10 -4 cts/s before then. 

From the data points that lie on the model curves, the 
simulated data, we synthesise the final data points by ran- 
domly perturbing the simulated data within the Gaussian 
errors. We simulated both single and broken power law light 
curves from all possible combinations of the parameters de- 
fined in section |2~T1 Samples of two synthesised XRT (SXRT) 
light curves are shown in Figures [T] & [2] 



2.3 Monte Carlo analysis: Fitting of data and the 
F-test 

The SXRT data, defined by given set of input parameters, 
are fit au tomatically using the simulated annealing method 
(§10.9 of iPress et alj|l992l and references therein) to min- 
imise the x 2 °f th e residuals. The SXRT data are first fit 
to a single power law, then a smoothly broken power law 
of fixed smoothing parameter, s, and the two are compared 
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t(s) 



Figure 1. A sample synthesised XRT (SXRT) data set includ- 
ing associated errors, overlaid on the underlying smoothly broken 
power law with a break at 1 day, smoothing factor, 5 = 5 and 
p = 2.4 (l>x > Vc,m)- 
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Figure 2. A sample synthesised XRT (SXRT) data set including 
associated errors, overlaid on the underlying single power law with 
p = 2.4 (i/x > I/ C ,m). 

via an F-test to evaluate their relative suitability. The data 
points are then re-perturbed from their original, on-model, 
simulated values and refit repeatedly (1000 trials) in a Monte 
Carlo analysis to get average values and la Gaussian devi- 
ations of the best fit parameters and F-test probabilities. 
This process is then repeated for all possible single and bro- 
ken power laws allowed by the possible combinations of the 
parameters in section [2~T1 



3 RESULTS 
3.1 The F-test 

Before we discuss our results, let us first consider the inter- 
pretation of the F-test, which returns an F-test statistic that 
we convert to an F-test probability, F pro b. In this situation, 
the F-test is a measure of the probability that the decrease 
in x 2 associated with the addition of the two extra param- 
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Figure 3. The distribution of the log of F-test probabilities, 
log(F prob ) for the SXRT data set described in the caption of 
Figure[T](p = 2.4, vy > v c ,m, ibrcak = 1 day, smoothing, S = 5). 
The vertical lines show the average value and la distribution, 
log(F prob ) = -8.3 ± 2.1 (as in Table |2)l of the data while the 
curve is the Gaussian with these values, for comparison. 



eters of a broken power law, a.2 and t brC ak, is by chance or 
not. As is common practice in light curve analysis literature, 
the following, though arbitrary, are approximately true: 

• Fp r ob > 10~ 2 =>■ favours single power law 

• 10" 5 < Fp ro b < 10~ 2 => neither favoured 

• Fprob i$ 10~ 5 => favours broken power law 

We presume a broken power law only if the probability 
of a chance improvement of the x 2 is very low (F pro b % 
10 ), while a single power law is presumed if the probability 
of a chance improvement is very high (F pro b > 10 -2 ). For 
the intermediate cases with moderate probability (10 -5 < 
Fprob % 10 -2 ) of a chance improvement, a single power law 
is usually presumed as it is the simpler model. An F pro b = 
10~ 2 corresponds to a 3a detection of a break while F pro b = 
10~ 5 corresponds to slightly less than 5a certainty. We have 
modified the F-test somewhat to return probabilities greater 
than 1 when the x 2 of the broken power law is greater than 
the single power law. If the underlying power law is single, 
a broken power law may give a worse fit as it is defined to 
break to a steeper slope only (a?i < ai), not a shallower 
slope as may occur during the data synthesis. 

When fitting the data, we use a number of different 
smoothing factors (s = 1,2,5, 10) and find that the F-test 
probability does not have a strong dependency, if any, on 
the smoothing factor of the model used to fit the data, only 
on the smoothing factor used to generate the data, 5*. For 
this reason we present only our F-test results for fits with 
a smoothing factor s — 5. The values of log(F pro b) given in 
Tables [2] [3] & [4] show that a true, underlying broken power 
law could be mistaken (—5 < log(F pr ob) < —2, bold font) 
for a single power law in a number of cases. Due to the la 
distribution of the returned F-test probabilities for a given 
set of parameters, it is clear that even in cases with F pro b J$ 
10 , a proportion could be misidentified as single power 
laws (Figure [3} . 

For a given value of p, those data sets where the un- 
derlying power law has a quite smooth transition from the 
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asymptotic values (S ~ 1), are most likely to be mistaken 
for single power laws. If the break is at the edge of the time 
range, it is also more likely to be misidentified and this is 
especially true in the case of a break at later times, when 
there are few data points post-break. There are a number 
of conflicting dependencies on the value of p itself. Firstly, 
the difference between the two temporal decay indices, Aa 
(p > 2): 

• Aa= (2+p)/4 (i/x > i>c,m) 

• Aa = (1 + p)/4 (v m <vx< v c , Wind) 

• Aa = (3 +p)/4 (v m <ux< Vc, ISM) 

so that a larger value of p implies a more pronounced, 
easily detectable break. However a larger value of p also 
causes steeper decay indices and fewer SXRT data points 
post break, making the break more difficult to confirm. If 
the break occurs at an early enough time (e.g., tbrcak = 
0.35 days) so that there is enough data post-break, this effect 
is unlikely to be an issue. The values of Aa are also depen- 
dent on the regime and environment, and we also observe a 
strong dependency on these. It is clear from the results for a 
constant density medium with v m < vx < v c (Table Q that 
there are far fewer misidentifications than in the results for 
a wind- like medium with v m <vx<v^ (Table [3}. This may 
introduce a bias in the detected breaks of GRBs as those in a 
constant density medium (z/ m < vx < v c ) are most obvious 
and hence easier to identify. 

Looking at the results for the single power laws in Ta- 
bles [2j [3] & [4] it is clear that if the underlying power law is 
not broken then log(F pro b) ~ or F pro b ~ 1 and it is highly 
unlikely to be mistaken for a broken power law, even con- 
sidering the distribution of the results. Given that the cases 
where intermediate probabilities (—5 < log(F pro b) < —2) 
are found all have underlying broken power laws, if a bro- 
ken power law does offer an improved fit then the underlying 
power law is likely broken. 

The above however, obviously depends on the inten- 
sity of the burst and dimmer bursts will have less certain 
conclusions, while the brighter bursts are less likely to have 
misidentifications. By normalising all our SXRT light curves 
at 1 day, we have made them as comparable as possible. We 
have simulated a number of representative bursts with phys- 
ical parameters, which show that in the circumstances de- 
scribed above, a broken power law may be 'hidden', i.e., an 
F-test will not favour the broken power law over the sin- 
gle power law, though the underlying power law is indeed 
broken. 

3.2 Return of parameters 

Unlike in the case of the F-test probabilities, which have 
little or no dependency on the smoothing factor, s, used 
to fit the data, the returned values of temporal slopes, ai 
& a.2 and break time, ibrcak do show some, though only 
marginally significant, dependency on this. To demonstrate 
this effect we show the case which gives a median value of 
Aa (vx > Vc,m, P — 2.4) and use a standard break time 
of 1 day (Table [5]l , though this discussion also holds for the 
other break times in our synthesis. 

We see from the best fit parameters returned from a 
given SXRT data set, that data fit with smoother breaks 
(s = 1) give later tbrcak, shallower ai and steeper ai than 



Table 2. Log of F-test probabilities, log(F pro b) for SXRT data 
with given values of p and vx > ^c,m- The SXRT data are 
generated from a single power law, and from smoothly broken 
power laws with given values of S, and breaks at 0.35, 1 and 
2 days. The cases where the F-test favours neither fit model 
(—5 < log(F pro b) < — 2), so a single power law would normally 
be assumed, are in bold font. 



tbrcak 





p = 1.8 


p = 2.4 


p = 3.0 


no tbrcak 










-0.10 ± 0.42 


-0.10 ± 0.31 


-0.07 ± 0.34 


0.35 days 








5 = 1 


-5.3 ± 1.8 


-4.4 ± 1.8 


-10.0 ± 3.2 


5 = 2 


-7.0 ± 2.1 


-9.6 ± 2.6 


-14.7 ± 4.3 


5 = 5 


-11.1 ± 2.4 


-16.5 ± 3.2 


-27.1 ± 6.3 


lday 








5 = 1 


-4.5 ± 1.6 


-6.7 ± 1.9 


-7.6 ± 2.3 


5 = 2 


-8.4 ± 1.0 


-11.7 ± 2.4 


-13.2 ± 2.8 


5 = 5 


-7.9 ± 2.0 


-8.3 ± 2.1 


-16.4 ± 2.8 


2 days 








5 = 1 


-4.9 ± 1.6 


-3.9 ± 1.7 


-3.8 ± 1.7 


S = 2 


-7.1 ± 1.8 


-4.8 ± 1.7 


-3.7 ± 1.6 


5 = 5 


-7.3 ± 1.9 


-4.3 ± 1.6 


-2.7 ± 1.5 



Table 3. Same as Table |2] but for the case that f m < vx < v c in 
a wind-like environment. 



tbrcak 





p = 1.8 


p = 2.4 


p = 3.0 


no tbrcak 


-0.13 ± 0.36 


-0.09 ± 0.38 


-0.08 ± 0.37 


0.35 days 
S = 1 


-2.4 ± 1.4 


-2.0 ± 1.3 


-5.8 ± 2.4 


5 = 2 


-3.4 ± 1.6 


-5.2 ± 2.0 


-7.7 ± 3.1 


5 = 5 


-6.5 ± 2.1 


-10.4 ± 2.7 


-17.4 ± 5.6 


lday 
5 = 1 


-2.3 ± 1.3 


-3.9 ± 1.8 


-4.3 ± 2.0 


S = 2 


-5.4 ± 1.9 


-3.6 ± 1.7 


-9.2 ± 2.8 


S = 5 


-5.4 ± 1.8 


-5.8 ± 2.0 


-13.4 ± 3.1 


2 days 
5 = 1 


-1.8 ± 1.2 


-2.5 ± 1.4 


-2.4 ± 1.6 


S = 2 


-3.1 ± 1.5 


-3.9 ± 1.7 


-3.4 ± 1.6 


5 = 5 


-3.8 ± 1.6 


-4.2 ± 1.8 


-2.6 ± 1.5 



those fit with sharper breaks (s = 5). Fitting the data with 
a smoothing factor, s = S seems to give the most accurate 
values of the returned parameters, though in the case that 
the break is smooth (S — 1), the returned values of espe- 
cially a2 and tbrcak are not very accurate. As an example, 
we can see from the case of S = 5 and s = 1 that fitting 
with a smoothing parameter less than the real smoothness 
(s < S) will underestimate ai while overestimating ai and 
tbrcak- The reverse is also true, in that if the data are fit 
with a smoothing parameter greater than the real smooth- 
ness (s > S, e.g., 5 = 1 and s = 5), ai is overestimated 
while ai and tbreak are underestimated somewhat. 

Whether or not s = S, there is quite a spread in the 
returned parameters for a given set of input parameters. 



Are the missing X-ray breaks in GRB afterglow light curves merely hidden? 5 



Table 4. Same as Table [2] but for the case that i/ m < ux < fc hi 
an ISM-like (constant density) environment. 



^brcak 





p = 1.8 


p = 2.4 


p = 3.0 


llO tl-irr.nl; 










-0.05 ± 0.56 


-0.08 ± 0.39 


-0.09 ± 0.34 


0.35 days 








5 = 1 


-8.8 ± 2.1 


-7.8 ± 2.2 


-19.6 ± 3.7 


5 = 2 


-10.5 ± 2.2 


-14.5 ± 2.8 


-22.4 ± 5.1 


5 = 5 


-15.1 ± 2.6 


-21.4 ± 3.2 


-35.8 ± 6.7 


lday 








5 = 1 


-6.7 ± 1.8 


-9.3 ± 2.1 


-11.9 ± 2.5 


5 = 2 


-10.2 ± 1.9 


-13.9 ± 2.3 


-17.2 ± 2.7 


5 = 5 


-9.2 ± 1.9 


-9.8 ± 2.1 


-7.8 ± 2.1 


2 days 








5 = 1 


-8.4 ± 1.7 


-8.9 ± 2.0 


-4.5 ± 1.9 


5 = 2 


-10.1 ± 1.9 


-10.1 ± 2.0 


-3.7 ± 1.7 


5 = 5 


-10.4 ± 1.8 


-9.3 ± 1.9 


-11.9 ± 2.3 



Table 5. The best fit parameters for a set of SXRT data for 
various values of fit smoothing parameter, s. The SXRT data is 
generated from a broken power law with p = 2.4 (fx > ^cm), a 
break at 1 day, and given values of smoothing parameter, 5. 







Ql 


a 2 


*brcak X 10 4 S 


input 


1.30 


2.4 


8.64 


5 = 


1 








s = 


1 


1.30 ± 0.18 


2.53 ± 0.37 


12.5 ± 8.3 


s = 


2 


1.38 ± 0.15 


2.32 ± 0.26 


9.5 ± 6.1 


s = 


5 


1.45 ± 0.11 


2.20 ± 0.28 


8.1 ± 5.4 


5 = 


2 








s = 


1 


1.19 ± 0.16 


2.66 ± 0.28 


11.6 ± 5.9 


.s = 


2 


1.28 ± 0.11 


2.44 ± 0.20 


9.3 ± 3.7 


s = 


5 


1.35 ± 0.08 


2.32 ± 0.20 


8.3 ± 3.2 


5 = 


5 








s = 


1 


1.20 ± 0.12 


3.27 ± 0.63 


14.6 ± 4.7 


s = 


2 


1.25 ± 0.09 


2.78 ± 0.40 


11.0 ± 3.1 


s = 


5 


1.29 ± 0.07 


2.50 ± 0.35 


9.0 ± 2.5 



This is similar to the result of Ijohannesson et al.l {2006), 
which suggests caution in the interpretation of results of 
broken power law fits to light curves. It could be an influence 
on the observed deviation of post-break slopes from a?2 = P- 



4 DISCUSSION 

We have shown that even an ideal, non-flaring, well sam- 
pled X-ray light curve synthesised from an underlying bro- 
ken power law can, in some cases, be as satisfactorily fit by 
a single power-law decay as a broken power-law. This would 
lead us, in the absence of further evidence, to assume a sin- 
gle power law and rule out the possibility of a break. This 
effect is most pronounced in bursts with late break times or 
high levels of smoothing ( S ~ 1), as many previously stud- 
ied jet breaks display fe.g. JZeh et al.ll2006l ). In the cases of 
high levels of smoothing, the breaks do not display sharp 
changes in the temporal decay index but a shallow roll-over 



from asymptotic values, well described by a smoothly bro- 
ken power-law, which makes accurate determination of the 
temporal slopes and break time difficult. 

The prototypical example of such a break is 
GRB 990510 for which well sampled B, V, R and I band 
light curves display an achromatic break (e.g.. lStanek et all 
1 19991 ) . This is accepted as a jet break even though the X- 
ray l ight curve as measured by BeppoSAX (Kuulke rs et al.l 
2000) is satisfactorily described by a single power-law, 
though the possibility of a b reak is not eliminate d. Likewise, 
in the case of GRB 060206 (jCurran et al.1 120071 ). the X-ray 
break is not well pronounced but supported by spectral and 
optical data. 

This has significant implications for the analysis of the 
numerous X-ray light curves that the Swift satellite has 
afforded us. For those X-ray light curves extending up to 
~ 1 day or longer, for which we do not have well sampled 
optical light curves, caution is required when making claims 
about the absence of breaks based solely on a comparison 
of the nominal fitted slopes. Multi-wavelength (optical & X- 
ray) temporal and spectral data up to late times are required 
to make a confident statement on the absence or presence of 
an achromatic break, and to determine the break parame- 
ters. This is particularly important when performing statis- 
tical analyses on a large sample of bursts, for making colli- 
mation corrected energy estimates, and for using GRBs as 
standard candles. 

Judging from the cases of GRB 990510 and 
GRB 060206, there may be a tendency, if not yet strongly 
significant, for X-ray light curves to have less pronounced 
or smoother breaks than optical light curves and more 
detailed theoretical models of jet breaks are necessary to 
clarify whether jet breaks could vary somewhat between 
wavebands. A full population synthesis, which is beyond the 
scope of this paper, is necessary to estimate the percentage 
of observed bursts that may have hidden X-ray breaks. 



5 CONCLUSION 

As underlying smoothly broken power laws may be hidden 
as single power laws, in XRT light curves, we should exercise 
caution in ruling out breaks based solely on a comparison 
of the nominal fitted slopes. We hence need to be cautious 
in implying chromatic breaks from optical and X-ray light 
curves where an X-ray break is not ruled out with a high 
degree of certainty. Multi-wavelength temporal and spectral 
data are required to make a confident statement on the ab- 
sence or presence of an achromatic break. There may be a 
bias towards detecting breaks from bursts with a constant 
density circumburst medium, as these are most obvious and 
easily detectable. 

We have shown that the fitted temporal slopes of 
smoothly broken power laws show significant variation about 
central values and that, especially in the case of a smooth 
break, the underlying values are difficult to extract via a 
fit. More accurate and wavelength specific descriptions of 
breaks, likely involving numerical simulations of the jet dy- 
namics, are necessary to better understand the observable 
signature of breaks and to clarify whether the breaks could 
vary somewhat between wavebands. 
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